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Abstract 

We study the self-energy amplitudes for the p — uj meson system produced 
by an effective field theory model at the quark level characterized by a finite 
range quark-quark interaction and an isospin breaking bare quark mass dif- 
ference rriu — md. Dynamical effects associated with confinement and dressing 
of quarks and the finite size of the produced qq meson modes are treated. 
Parameters of this approach, previously constrained by soft pion physics, are 
here supplemented in a minimal way to reproduce and gfym^. Calculations 
are carried out up to the one pion loop level and the predictions for 
mixing, mass splitting, the p width and the symmetry breaking wvrvr cou- 
pling constant and form factor are presented. The contribution of the direct 
oj TTTT process relative to the a; vrvr process is discussed within this 

model. 

PACS numbers: 24.85. +p, 12.40.Yx, 12.39.Ki, 13.75.Lb 
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I. INTRODUCTION 



Recent work on the pion and kaon systems has shown [|T],3 that a very efficient description 
of masses, decay constants and electromagnetic form factors in terms of QCD degrees of 
freedom can be implemented through a parameterized form of the dressed quark propagator 
used in conjunction with truncations of the coupled Dyson- Schwinger equations (DSE). 
Chiral symmetry considerations for these systems eliminate the need for explicit information 
about gluon propagators and the gluon-quark vertex apart from that which is implicitly 
provided by the dressed quark propagator. At the simplest level this DSE approach to the 
QCD modeling of hadron physics has been successful in a number of applications involving 
the Goldstone bosons including vrvr scattering and the anomalous 7*717 ■ynmi [|| and 
■jTip processes. The dynamical format of these initial studies can also be generated 
from a model field theory with QCD degrees of freedom containing an effective finite range 
gluon propagator. In particular, a treatment at the O*'^ order in the loop expansion of 
the produced qq Goldstone meson modes will yield the same dynamical structure with a 
consistently specified content for the various dressed elements. We take this latter point of 
view for the present work and consider a minimal extension beyond the Goldstone boson 
sector to address some interesting issues in the p — u vector meson system. A previous work 
on p — u mixing with this approach, |^ limited to 0*'' order in meson loops, is here extended 
to the one pion loop level to address a wider set of phenomena. 

Quark confinement should be an important element in the bound state dynamics and 
hadronic decays of the p — uj system. Otherwise, with typical constituent u and d quark 
masses of about 350 MeV, the nearby spurious qq production threshold would have undue 
influence. The approach we take incorporates quark confinement through the absence of a 
mass-shell pole in the employed dressed quark propagator. The necessary momentum de- 
pendence of the quark self-energy amplitudes implies a finite range for the effective gluon 
propagator. The associated quark momentum dependence of the meson-quark vertex func- 
tions reflects the finite size of the qq correlations. To guide our development of an effective 
action for vr, p and uj that preserves these realistic features we use the two flavor version of the 
Global Color Model (GCM) []TU|,|rT[. The employed dressed quark propagator model is con- 
strained by soft chiral physics and is a simplified version of what is currently available [jl|,|^ 
in order to facilitate the type of calculations necessary for this work. Information about the 
dressed gluon propagator, necessary to treat mesons that are not Goldstone mesons, is here 
introduced in a minimal way that is constrained by calculated values of m^^ and Qp-j^-n. 

The zero range limit of the model we employ is very closely related to the Nambu- 
Jona-Lasinio model (NJL) ||12| which has been extensively used to model hadron physics. 
However, in contrast with the present approach based upon the GCM, many important loop 
integrals encountered within the NJL model must be regulated via additional parameters. 
Also quark confinement cannot be generated from the contact NJL interaction and the 
dynamically generated quark mass is a constant (constituent mass) rather than the (running) 
mass function which is a necessary consequence of finite range gluon propagation. The 
associated nonlocalities have a profound effect upon the nature of the calculations and the 
role that is assigned to phenomenology. A previous study of the p — uj system within the 
NJL model was limited to the quark loop level for the mixing amplitude, but did generate 
the p propagator, including the 27r width, at the pion loop level [T^. In Refs. |ll^ , p!5[| an 
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extended NJL model that includes an extra confining interaction has been applied to the 
p — uj system in a way that is the most comparable to the approach we take. 

One of the topics we investigate is p° — mixing generated by the bare quark mass 
difference — rriu- Recent interest in a description oi — uo mixing from the quark level 
has centered around the momentum dependence of the mixed self-energy amplitude and the 
possible contribution to the charge symmetry breaking (CSB) component of the nucleon- 
nucleon force. |^,|TB|,|T^ The quark loop mechanism produces a momentum dependence that 
suppresses the amplitude at the small space-like momenta appropriate to the mixed p — uj 
exchange mechanism. Here we investigate the one pion loop correction and find that the 
space-like mixing still remains negligibly small. A more promising explanation for the CSB 
nuclear force has been found in the meson-nucleon vertex 



Our study of the pion loop contributions to the self-energies of the p — uj system produces 
calculations of the yovrvr and cuvrvr coupling constants and form factors as well as the 2% 



width of the p and the mass splitting 



m„ 



Initial studies within the GCM of 



T9[| and Qpt^t^ ||2^ approximated hadron mass-shell information by behavior at zero 



meson momenta. This was to avoid ambiguities associated with numerical extrapolations 
of dressed quark propagators to the complex momenta needed to put hadron momenta on- 
shell in this Euclidean space approach. The correct mass-shell definitions can be maintained 
in the present work due to the assured analytic properties of the model quark propagator 
employed. 

The amplitude for the direct process uj 
out in the analysis of e"*" e" 
amplitude |^ . 



TT vr is conventionally assumed to cancel 
vr vr data to extract a value for the physical p^ — uj mixing 
Here uji is the pure isospin component of the physical p. However, it has 

vr vr amplitude, the 



recently been argued that, without a theoretical model for the uji 
error on the extracted mixing amplitude is much greater than previously thought. Evidence 
for a significant contribution from the cj/tttt process is obtained in a recent QCD sum rule 
analysis of the mixed-isospin vector current correlator ||2^. In contrast, the value for Quij-k-k 
recently estimated from the extended NJL model of Ref. [|15| suggests a negligible effect. 
In the present work, the strength of this G-parity violating cj/tttt coupling, as driven by 

— in the underlying quark loop, is found to be significant enough to require an 
increase of typically 20% in the required mixing amplitude. 

In Sec. |T| we describe the underlying quark model action and also general aspects of the 
resulting effective meson action that we shall use. The connection between them, i.e. the 
generation of qq meson modes from the GCM action, has been covered in previous works 
and key elements are here reviewed briefly in the Appendix. Specific aspects of the pion 
sector, and the model dressed quark propagator which is constrained by it , are described in 
Sec. |T|. The mixed p — uj meson sector produced at tree level by integrating out the quarks 
is described in detail in Sec. At that same level, the produced pirir and ujtttt interactions 
are identified in Sec. |V] and the calculated coupling constants and form factors are described. 
The pions are integrated out up to the one loop level in Sec. and the contributions to the 
mixing as well as to the p width and mass shift are described. Also presented in Sec. [V| is 
a discussion of the relative contribution oi p^ — uj mixing and the direct process ujj ^ tt ir 
to the physical coupling strength g^jnn- Further discussion follows in Sec. |VII. 
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II. QUARK MODEL 



To obtain a practical covariant description of both meson substructure and interactions in 
terms of dynamically dressed quark degrees of freedom we use the Global Color Model |lTO| , p!T . 
This is an effective field theory at the quark level described by the action 

S[q, q] = I d'x g(a;)(7 ■ + m)q{x) + ^ | dVy r^{x)D,,{x - y)jM (1) 

where the quark currents, j^{x) = ^(a;)^7^g(x), interact via a phenomenological gluon two- 
point function {D^^{x) = 5^yD{x) in the Feynman-like gauge chosen here). We consider 
only two flavors {u and d) and use a Euclidean metric throughout, such that a ■ b = a^b^ 
and {7^, 7i.} = 25^,^, with 7^ = 7^. The bare masses are represented by the diagonal matrix 



m = {mu,md)- The closely related NJL model [0 is specified by the zero-range limit of 
Eq. (|I]) together with a set of coupling constants for the independent meson modes. 

There has been extensive work on both the meson and baryon sectors that arise from 
hadronization of the GCM and a comprehensive review is available. For just the lowest 
mass modes, the meson action for the effective local field variables vf, can be written 



as 



00 / \n 

S[tt, p, cj] = Tr ^ [5*0(275^^- Trf^ + i'y^ujfJv + ilf^r ■ pj^v)]"' 

n=2 ^ 

J ^ 2D{x-y) ^ ' 

For completeness, we collect the main steps in the transition from Eq. (|l|) to Eq. (0) into 
the Appendix. In Eq. (H) Tr denotes the trace over color, flavor and spin as well as the 
continuous space-time variable. The quark propagator is S'o(p) = [i '^AijP') + B{p^) + m]~^ . 
The first term of Eq. @) represents the sum of all possible meson couplings mediated by a 
single quark loop and the second term can be described as a contribution to the meson mass 
terms. 

For simplicity, we have made the common approximation in Eq. (0) that each meson 
Bethe-Salpeter (BS) amplitude is truncated to its single canonical Dirac matrix covariant 
with r being the associated scalar function. In studies of the BS equation that allow a more 
complete selection of covariants, p6| , p7| this is found to be the dominant component. In 



Eq. d^) we have also used the result that for the pure isospin components = F^ which fol- 
lows from the ladder BS equation when coupling terms generated from — are ignored. 
This defines the vector meson states of good isospin whose mixing we shall subsequently 
study. In a momentum representation the meson-quark coupling implicit in Eq. (Q) is, e.g., 
q{q^)i'^fj,ujfj_{P)Vv{q] P)q{q-) where the quark momenta are g± = g ± P/2 with P being the 
meson momentum. Thus Fy(g;P) plays the role of a form factor but has, in principle, a 
well-defined dynamical content defined by the model at the quark-gluon level. In this work 
we utilize the dynamical structure of Eq. (0) to combine previous work on the pion sector 
with a phenomenological approach to Fy to explore whether important features of the p — uj 
system can be described. 
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III. THE PION SECTOR 



The free pion part of the action in Eq. (Q) is 

1 r d^P 



where the inverse propagator for the effective local pion field is 

K\P') = tr / ^[^o(g-)^75r^o(g+)^75r] f ^(g; P) + ^ / (4) 

Here g± = g ± ^ and r is any single component of r. Throughout this work the symbol 
tr denotes a trace over color, flavor [Nf = 2) and spin. We note that is a functional 
of r,r and that 5A~^/5r^ = is the ladder BS equation (projected onto the i'j^T channel) 
for the amplitude T.„. To produce r,r that way requires a specification of the model gluon 
propagator. However an indirect route is provided through chiral symmetry by the Goldstone 
theorem in terms of the dressed quark propagator. The realization relevant to the present 
treatment is that, in the chiral limit, the ladder BS equation for the amplitude r^(g;P) 
becomes identical to the ladder Dyson-Schwinger equation for B^^q^) = B{q'^;m = 0), the 



scalar part of the dynamically generated quark self-energy |28[. Thus the chiral limit pion, in 



this approach, is automatically both a massless Goldstone boson and a finite size qq bound 



state. For a small bare mass, numerical studies reveal that is still well approximated 
by the scalar quark self-energy amplitude. 

In this work we take Tniq; P) ~ B[q^] m) and use mav = {rriu + trid)/'! as the bare mass 
for both the u and d. Near the mass shell, we have the general form A~^(P^) = (P^ -|- m^)/^, 
where is the pion wavefunction renormalization constant. The (constant) second term of 
Eq. @) may therefore be eliminated if we require that the GMOR relation in the form 
A^^(O) = 2mav\{qq)\ be satisfied. Here {qq) is the isospin average vacuum condensate. In 



this regard we note that recent investigations [^,0] find that a somewhat differently defined 
quantity is better suited to this relation. The numerical consequences for the present work 
are not significant. Our procedure yields 

k-\P^) = 2m,.|(gg)| - tr j ^[5o(g_)5o(-g+) - So{q)So{-q)] B^q'), (5) 

where we take \ {qq) \ = ^ tr '^^^oiQ) with the scale specified by the cutoff //^ = 1 GeV^. 

Now is fully specified in terms of the quark propagator. Numerical evaluation of Eq. (^ 
identifies the pion mass through the zero at P^ = —Tn\. 

We also identify Z^{P'^) defined by A-i(p2) = (p2 ^ ml)Z^{P'^). The momentum de- 
pendence of Zt.^P'^) implies a dynamical mass function mT^^P"^) due to qq substructure. We 
calculate from = Z^^^P^ = —m^). A physical normalization (unit residue at the mass- 
shell pole) is produced by absorbing into the pion fields. For convenience, we redefine 



the fields to absorb ^J Zj.i^P'^) so that the propagator becomes A^^ = (P^ + m^) ^. The 
associated vertex amplitude is then the dimensionless quantity F^(p; P) = B{p'^)/ J Zt,{P'^). 
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At the mass-shell this procedure satisfies the canonical Bethe-Salpeter normalization con- 



dition |32[. For any pion momentum, this choice defines the effective pion -quark vertex 
and related pion field consistent with the point-pion propagator = (P^ -|- m"^)'^ to allow 
contact with common procedures in nuclear physics. In various applications we test the role 
of the composite nature of the pion in by replacing the function Zt^IP'^) by its mass-shell 
value m^). Such tests will be described later. 

For the dressed quark propagator, we use the representation So{p) = —i'j ■ p (7v{p'^) + 
(Js{p^)i and employ the following simplified form of previously developed parameterized 
amplitudes n 



as{x) = c(m) e-^^ + - (l - e'^^) , (6) 



2x — 1 + e~^^ 

^v{.x) = — c(m) m e"^'', (7) 

with X = p'^/X^, cTs = XcTs, o'v = A^cry, ^ = rn/\, where m is the bare quark mass 
and A is the momentum scale. This So{p) is an entire function in the complex momentum 
plane, a condition sufficient to ensure the absence of quark production thresholds in S- 
matrix amplitudes for physical processes ||3^. The original parameterization is guided by 
a confining model DSE and by the behavior found in realistic DSE studies [@,^. It 
is consistent with pQCD in the deep Euclidean region apart from In(p^) corrections. The 
parameters are A = 0.889 GeV, c^rfiav) = Cq = 0.581, niav = 16 MeV. The soft chiral 
physics quantities produced by these parameters are = 90.1 MeV, = 143 MeV, and 
(qq) = {173 MeV)^. The more elaborate original parameterization ||l[ provides a better fit 
to the above and also an excellent description of the pion electromagnetic form factor and 
the TTTT scattering lengths. 

With the isoscalar quark propagator set as above, the isovector term is generated from 
the bare mass dependence by substitution of m = rhav + y^m with 6m = rhu — frid- The 
constant c(m) mainly represents the strength of dynamical chiral symmetry breaking that 
generates the dynamical mass function, M(p^) = crs{p'^)/<^v{p'^)- However a weak bare mass 
dependence for c{rh) is necessary to simulate the bare mass dependence of realistic DSE 
solutions. This leads to a small isovector component of c(m) and we write c(m) = cq + T3C1, 
where cq = c^friav) is given above, and ci = (c^ — Cd)/2. To determine the latter we fit to 



the bare mass dependence of as{p = 0) produced by realistic DSE studies Only a 

linear behavior ci = d6m/2 characterized by the slope d = c'{rh) < is important, and with 
Sm = —4 MeV, we have Ci = 0.0048. This 1% isospin-breaking effect Ci/cq is of the order 
of magnitude expected. The quark propagator can be written as So{p) = Sq{p) + t^Sq^p) 
where the latter isovector component is responsible for the various — u mixing mechanisms 
considered here. 



IV. THE p-u SECTOR AT TREE LEVEL 

We treat the transverse modes in the p and uj channels and use fields that contain the 
transverse projector T^u{P) = V ~ PimPu/P"^- That is, p^,{P) = T^^{P)py{P) and uj^{P) = 
Tfj,u{P)uJu{P) . In a matrix notation where denotes (p^,^^^), the tree-level action from 
Eq. (H), up to second order in the fields, can be written as 
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(27r)^ l 



(8) 



Here, in meson channel space, is diagonal and the only non-zero elements of W are the 
off-diagonal ones that provide — uo coupling. The superscript q identifies the quark-loop 
contribution to distinguish from a different contribution introduced later. The diagonal 
inverse propagator is given by 



A-!(P) 



tr 



tl{q- P) [So{q-)il,fS,{q+)i^J] + 9 5, 



D r) 



(9) 



Here / are the relevant flavor factors (1/ for the u, and r for the p) and q± 



g ± ^, where 

P is the meson momentum and q is the loop momentum, or equivalently the qq relative 
momentum. The amplitude Ty is defined in the diagonal sector by the ladder BS equation 
produced by SA'^/STy = 0. Thus Ty contains no isospin mixing mechanism and is the 
same for p and u. The diagonal p and u elements of the inverse propagator are likewise 
identical. The transverse component A~^(P^), defined by A~J(P) = T^,^(P)A~^(P^), is 
given explicitly by 



A-\P' 



24 



+ 9 



(2^ 
D r) 



p2 2 



2(P ■ qf 
3P2 



(10) 



Here we have ignored a small contribution from the isovector component of the quark prop- 
agator by using niav as the bare quark mass. 

The mixed self-energy, H^^(P), is produced from the same quark loop term of the effective 

action of Eq. (g) that produces A^J. The mechanism is illustrated in Fig. 1. The explicit 
expression is 



dS 
(2^ 



[So{q-)hi,r3So{q+)hu] f y(g; P), 



(11) 



which is identical to the first term of Eq. (^ except that one of the flavor factors is Ts and the 
other is If. Thus, n^^,(P) is nonzero only if 5*0 has an isovector component. The structure 
of Eq. ([TT|) is corresponds to first-order perturbation theory for isospin breaking effects in 
that the internal meson states or BS amplitudes that enter are from the unperturbed or 
isospin symmetric sector. The mixed self-energy may be expressed as the difference of quark 
loops and the leading contribution is proportional to Sm. The transverse component is 
fl«(p2) = fl„ - fid, where 



n/(P^ 
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d^q 
(27" 



P)\aUql)4{Q'- 



P^ q^ 2(P ■ qf 
4^3 3P2 



(12) 



To summarize the transverse p^ — uj sector, the inverse propagator obtained at the quark 
loop level may be written as 
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^ fi.(p2) A-i(p2) j • ^^^^ 

We shall refer to the basis here as (p?, oji) to indicate pure isospin states. 

The composite qq nature of the mesons is reflected in the fact that the diagonal elements 
contain a dynamical mass function. The general form 

A-i(p2) = {P^ + ml)Z^{P^) = (P^ + mUP^)) (14) 

identifies the degenerate mass my of the pj and ui states. It is convenient to carry out 
part of the necessary renormalization at this stage. Without the off-diagonal terms of 
Eq. (p!3D, a physical normalization (unit residue at the mass-shell pole) would be produced 
by absorbing at least the on-mass-shell value of ^/Z^ into the fields. As with the pion sector, 



we choose to absorb the function yZi(P^) into the fields so that the resulting unmixed 
propagators have the standard point meson form, to The associated BS amplitudes are now 



Tviq',P) = Tv{q', P)/ J Zi{P'^), which are dimensionless quantities that, at the mass-shell. 



satisfy the standard normalization condition [^] for a BS amplitude. We note that the 
meson-gg coupling strength (magnitude of Fy at the meson mass-shell) is thus a prediction 
of this approach. The quark loop contribution to the — u sector is now 

where n'^(P^) = II^(P^)/Zi(P^). The role played by the dynamical qq substructure of propa- 
gators may be obtained by holding the function Zi(P^) at its on-mass-shell value to simulate 
a structureless meson. 

To facilitate these exploratory calculations, we adopt the simple form ry(g;P) = 
exp{—q^/a?) which has previously been found to provide an adequate one-parameter 
representation of the relevant BS amplitude in a variational approach [^]. The physical 
amplitude Fy is independent of N . The choice of effective gluon propagator D{r) in Eq. (|10|) 
is made indirectly by ensuring that the resulting constant value of the second term produces 
a zero in A~^(P^) at P^ = —my where my = = 782 MeV. The range parameter a is 
chosen to reproduce the empirical value of gpT^T, as discussed in the next Section. 

The transverse mixing amplitude 11^ (P^) obtained from this approach has been discussed 
in a previous work ||^. The present parameterization of the input to this model is slightly 
different but only small quantitative changes occur in the calculated 11'^ (P^) and the con- 
clusions of our previous work remain unchanged. 



V. THE pvrvr AND ojtttt INTERACTIONS 

From the effective action in Eq. (|^), we identify the interaction term 

5' [ptttt] = -Tr Soij^T ■ p^TviSoh^T ■ nT^f . (16) 

Use of the momenta shown in the quark triangle diagram of Fig. ^ yields the explicit ex- 
pression 
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S [pnn] = z / p,{Q) . n{-P - |) x vf(P - |) A^(P, Q), (17) 



where the vertex is 

A^(^, Q) = J 0^,^v(k + ^; Q)r.(fc + 1; - |)r.(A: - f ; ^ - ^> (is) 

with the discrete loop trace denoted by 



T^{k,P,Q) = ti 



Soik + ^ + ^) i-ff,Soik + ^ - y) h5So{k - ^) ij^ 



(19) 



The symmetry properties of A^(P, Q) that follow are A;^(P, Q) = -A^(-P, Q) = A;^(P, -Q). 
Hence the most general form is 

A^(P, Q) = -P^ Fp^^{P\ Q\ (P ■ Qf) -Q,P-Q Hp^^{P\ Q\ (P • Qf). (20) 

With both pions on the mass-shell, {P — = (P+fi)2 = — m^. Equivalently, P-Q = 
and P2 = -ml - ^ so that the simplified form of Eq. (^) is A^(P,g) = -P^Fp^^{Q'^). 
The coupling constant is defined as the mass-shell value (7^,^^ = Fp^^{Q'^ = —m?p). If the 
form factor is held at this value for all momenta, one obtains the point coupling limit 

5 [ptttt] = -igp^^ I P, p,{Q) . 7r(-P - |) x n{P - |). (21) 

which is equivalent to the more standard form 

5* [pirn] = -gpTTTT j d'^x pp{x) ■ 7r(x) x d^TT{x). (22) 

With the parameterized BS amplitude ry(g;P) oc exp{—q^/a^), we take a = 0.194 GeV^ 
to reproduce the value g'^^^ = 6.05 as inferred from the experimental p decay width 

of 151 MeV. The calculated form factor FpT^T^{Q'^) is displayed in Fig. ^ for timelike and 
spacelike momenta in the vicinity of the mass-shell. A previous study |^ within the GCM 
made the approximation gpT.^, ^ Pp^^(P = Q = 0) to avoid the occurrence of complex 
momenta in the arguments of propagators and vertex functions in the quark loop. From 
Fig. ^ it is evident this approximation can underestimate the value by almost a factor of 
two. 

The pion loop mechanism to be considered later requires the ptttt vertex for off-shell pion 
momenta. Only the form factor Fp^^^^ is relevant due to the transverse p condition. Results 
from numerical evaluation of Eq. (|T8|) over the relevant range of P^ and for P ■ Q = 
are displayed in Fig. ^. For later use we parameterize the results in the form 

i^p..(P^Q')^/p..(g')e-^^/^?('5^^ (23) 

The constraint P ■ Q = is applied to facilitate subsequent use in the tt loop integral and 
it is exact at the important on- mass-shell point for both pions. The functions fp-K^iQ'^) 
and Ap(Q^) are chosen to provide a fit to the calculated p momentum dependence of the 
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vertex to within 1% for between — and 0.6 GeV^. We use fprnriQ"^) = J2n=o'^nZ^ and 
'^p(Q^) = J2n=o^nZ"' with z = Q^/yU^ and = I GeV. The parameters are (00,01,02) = 
(3.04, -5.20, 2.37) and (60, • • • , ^e) = (1-28, 0.252, 0.110, 0.171, 1.08, 1.80, 0.865). 
The tuvrvr interaction term from the effective action of Eq. (|^) is 

S[ujnn] = -Tr [S'oi7^t^^ry(S'oi75r ■ Trr^)^] . (24) 

This is identical in form to the pvrvr interaction except for the isospin structure. If the 
only isovector quantities involved here are the pion fields, the resulting expression is just 
Eq. ([T7|) but with isospin structure vf ■ vr. The symmetry properties of AJ^(P, Q) then show 
that the result is identically zero as required by G-parity conservation. With extra isovector 
components provided by the quark propagators, the G-parity violating result is 

SsIootttt] = I J u,{Q) 3 • 7f(-P - Q/2) x 7r(P - Q/2) K{P- Q), (25) 

where the vertex is 

K^P. Q) = j 7Sirv(fc + ^; g)r.(fc + ^:-P- |)r.(A: - f ; - ^)N,{k, p, q). (26) 



(27r) 

Here, to lowest order in isospin breaking, the loop trace is proportional to 5m and is given 

by N^, = - A/'OO^ + A/'Oi°, where 



iVf^(fc,P,Q) = tr 



(27) 



with the isospin components of the propagators defined by = + ts^q. 

The symmetry properties of A||^(P, Q) are identical to those for hPjyP-i Q) and hence the 
general form is 

A;i(P, Q) = -P^ F^^^{P\ Q\ (P • Qf) -Q,P-Q H^^P^ Q\ (P ■ Qf). (28) 

At the triple on-mass-shell point, numerical evaluation of PajTrw yields the coupling constant 
Qujj-K-K = 0.105. This is almost five times larger than the value 0.0236 recently obtained in 
Ref. [|r§]. The 2tt width of the physical u is 0.185 MeV which implies g^^* = 0.21. The 
theoretical quantity Quij-k-k obtained so far refers to the pure isospin uj state and a small 
admixture of pj can make a significant contribution to the physical u ^ 2n decay. We will 
return to this topic later. The calculated off-shell form factor F^t^t,{P^,Q^) is displayed in 
Fig. 1^ for P ■ Q = 0. A common assumption is that the pirn and ujtttt form factors simply 
scale by the ratio of their coupling constants so that the pion loop p^ — uj self-energy may 
be estimated by scaling the corresponding p self-energy ||2T1J36| . From Figs. § and ^ the 
present calculations show that this can be qualitatively reliable for momenta where the form 
factors are significant. The departures that do occur are due to the different momentum 
dependence of the isoscalar and isovector components of the quark propagator. 

For the subsequent pion loop calculation we parameterize the ujtxti form factor in the 
same way as in Eq. (^. That is, fujnniQ'^) = En=o'^n^" and A^(Q^) = p En=o^n^'' with 
z = Q^/p"^ and p = I GeV. The parameters are (oq, Oi, 02) = (7.07, —6.91, 3.21) x 10^^ and 
(60, ...,be) = (118, 6.76, -1.72, -0.391, -0.172, 0.109, 0.140) x lO^^. 
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VI. THE p-LO SECTOR UP TO ONE-PION-LOOP 

The relevant low-order terms in the effective action for vr, p, uj that we have produced by 
integrating out the quarks are 

S [vr, p, uj] = S2 [p, iv] + S2 [vr] + S[p7ni] + SIluttti] + ■■■. (29) 

Here the second-order terms 5*2 [p, uj] and §2 [tt] represent respectively the p — u sector con- 
taining the quark-loop mixing mechanism and the free pion sector. To integrate out the pion 
fields, it is convenient to first combine the last three terms of Eq. ( pQ] ) so that it becomes 

S[7r,p,u] = S[p,u] + y ^^n,{P')Dr,\P',P)7r,{P). (30) 

The three terms in 

DT^\P', P) = A-^iP', P) + V,,{P\ P) + W,,{P', P) (31) 
correspond to the free pion term given by 

A-.i(P', P) = (2vr)%5^(P' + P)A-\P% (32) 
the pTTTT interaction term given by 

V.,iP',P) = 2ie,,,pl{-P' - P) Ki-^^-^ -P' - (33) 
and the uinr interaction term given by 

W,,{P',P) = 2ze,,3^,(-P' - P) _p' _ p). (34) 

Integration over the pion fields allows a new effective action S[p, uo] to be identified from 

Z = N J DpDuDn exp [-S[7r, p, cu]) = N' J DpDu exp (-5[p, to]) (35) 
by using the functional integral result 

/ Dnexp (^-^ / ^j^MP')DT^\P',P)7r,{P)^ = exp (-^TrLn D'^^ . (36) 

We absorb the field-independent term ea;p(— |TrLn A^^) into the normalization constant to 
arrive at 

S [p, Lo] = S2 [p, u] + ^TrLn(l + A{V + W)). (37) 

The second term here gives the vector meson coupling to all orders produced by a single 
pion loop. 

11 



The quadratic part of the second term of Eq. (|37D adds a contribution to the vector 
meson inverse propagator to that given previously by Eq. (|15D at the quark loop level. In 
the transverse — uo sector, the net result for the inverse propagator is 

^ ^^>-\ n^-(p2) P^ + ml + W^{P^) )' 

where the net mixing amplitude is n'"^(P^) = n'^(P^) + n'^(P^). The additional terms pro- 
duced by the pion loop are the diagonal elements I\.pp{P'^) and n'^'^(P^), and the mixing 
term 11'^ (P^). From Eq. ([37|), the self-energy contribution I\.pp{P'^) is given by the transverse 
component of 

n^^.(P) = -4 / ^A.(g+) A.(g_)A^(g; -P)K{<1: P), (39) 

where q± = q ± P/2, and the summation over pion isospin labels has been carried out. This 
contribution is the same for each isospin component of the p. Similarly, the mixed p^ — uj 
self-energy contribution generated by the pion loop is given by 

K^{P) = -4 / ^A.(g+) A.(g_)A^(g; -P)KW. P). (40) 

and n'^(P^) is the transverse component of that. The uj self-energy contribution is obtained 
from Eq. ( p9|) by the replacement of by A'^. It is quadratic in the small symmetry 
breaking mechanism and is ignored here. 

In the absence of mixing, the real part of Upp^P"^) generates a mass shift for the isospin 
eigenstate p/. Since the corresponding quantity in the uji channel is much smaller, and 
neglected in this work, this mass shift represents the p — uj mass splitting. For timelike 
momenta such that P^ < —4m'^, both n''''(P^) and n'^(P^) have imaginary parts associated 
with the decay p — > 2tt. The imaginary part of n''''(P^) produces the 2n width of the pi 
state. 

These observations are, in principle, modified by the presence of nonzero off-diagonal ele- 
ments in Eq. (|38|) . To see this, we construct the eigenvalues of the inverse meson propagator, 
P~^(P^). These are, to lowest order in the total mixing amplitude, given by the expres- 
sions Ap(p2) = p^ + ml. + WP{P^){1 + e2(p2)), and X^{P^) = P^ + ml - WP{P'')e''{P^)), 
where 

e(p2) ^ -iL-kLi (41) 
I{PP{P^) 

is a natural measure of off-diagonal coupling. Since the typical size is |e| ~ 0.03 in the 
mass-shell region, we shall ignore the second order mixing effects on the mass shifts. In the 
diagonal basis we therefore have, through first order in mixing, 

p. :,„.). (42) 

where we have made the identification = my. We identify the mass of the (mixed) p 
eigenstate from the position of the zero in the real part of the upper eigenvalue in Eq. (^2]). 
That is, 



12 



(43) 



With the definition f(p2) 
be written as 



-lmIl''''{P'^) /nip, the inverse propagator in the p channel can 



V~\P') = (P2 + ml- impTiP^)) Z2{P^ 



(44) 



where r(P^) = f (P^)/Z2(P^). The function Z2(P^) arises from the momentum dependence 
of the real part of the pion loop self-energy. The on-mass-shell value, which is to be absorbed 
into the p field as a renormalization constant, is given by 



1 + Ren^^ i-ml 



(45) 

with the prime superscript denoting differentiation with respect to the argument. We choose 



to absorb the function ^ Z2{P^) into the field P/i(P) so that the resulting propagator has 
the conventional Breit-Wigner form. The physical width is given by 



Imn''''( 



mpZ2{ 



(46) 



Our numerical results for the p mass shift and width are obtained from Eq. ( ^9|) in the 
specific form 



(g2 _ (g . Pf/P^)e 



(47) 



where the parameterized pim vertex functions from Eq. ( p3D have been used. After reduc- 
tion to a two-dimensional integral it is convenient to employ the Feynman combination of 
denominators to obtain 



4 



da / dss 



+1 



dzVl 



-2{s-2azsP+a^P'^)/Xl{P'^) 



2 • 



(48) 



The numerical results for Re n^^(P^) and Im TI.pp{P'^) are shown in Fig. ^. The self-consistent 
solution for uip produced by Eq. (|^) is displayed in Fig. ^ in the following way. The 
quantity — (P^ + m^) is plotted as a long dash line and its intercept with the solid line 
representing Ren''''(P^) identifies the mass shell point. We use = my = 782 MeV and 



find nip = 761 MeV giving 



m„ 



-21 MeV. The experimental value is —12.0 ± 0.8 



MeV. At this mass-shell point we find UPP{-mf,) = -0.031 - O.lli GeV^ The calculated 
field renormalization constant from Eq. (^) is Z2(— m^) = 0.90. Then Eq. ( ^61) yields the 
width Tp = 156 MeV while the experimental value is 151 MeV. 

Fig. 1^ also illustrates the contribution to the mass splitting result made by the momentum 
dependence of the meson self-energies generated by the qq substructure. This is most marked 
for the pion and enters the pion loop integral in Eq. (^) through the factors Z~^(g^) which 
are contributed by each pion propagator. (In our calculational procedure these factors have 
been moved into the vertex functions for convenience.) If each function Z^^ is held constant 
at its mass-shell value (as would be the case for a structureless pion), the full calculation 
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in Fig. (solid line) becomes the dot-dashed line. The intercept then indicates that the 
mass shift would be essentially zero. In contrast to this, the short dashed line indicates that 
the influence of the corresponding quantity Zi(P^), from the substructure of the p, is quite 
negligible here. This is because the range of variation in is very small. 

In a parallel manner, the pion loop contribution to the transverse p^—u mixing amplitude 
is calculated from 



w(p^^- rp2w ,p2. r A (g^-(g-p)VP>'^^-^- > 

6 J {111} (^q jf. Jiy rnj^ - (^g - rnj^ ^ le 



where A ^ = + A^^^. The results are shown in Fig. |[ An experimental value ||22 



for the total real mixing amplitude is also displayed. The pion loop mechanism produces 
a significantly smaller real amplitude than does the quark-loop mechanism and does not 
have a node near zero momentum. In Fig. ^ we illustrate the contributions to Re 11'^ made 
by the dynamical nature of the meson self-energy functions. When the pion quantities 
Zj.i^q'^) are held fixed at the mass-shell value, the full calculation (solid line) reduces to the 
dot-dashed line. As with the p self-energy, the dynamical qq substructure of the pion has a 
considerable influence on the pion loop mixing mechanism near the vector meson mass-shell. 
The difference between the solid and dashed lines shows that the corresponding vector meson 
substructure plays a lesser role. We compare the quark loop (dashed line) and pion loop 
(dotted line) contributions to the real mixing amplitude in Fig. 0. The solid line is the sum 
of the two. At the mass-shell, the quark loop dominates with the pion-loop contributing 
about 20% of the total. For spacelike momenta, the pion loop contribution is negligible. 

The change of basis associated with the diagonal form of the inverse propagator given in 
Eq. ( ^21) produces physical fields given by 

(50) 

Currents that couple to pi and uj are mixed in the same way. Thus the resonant p and lo 
contribution to the time-like pion charge form factor is given by 

F^^iP') = F^UP')V^iP')U,iP') + F,UP')V,iP')f,,iP'), (51) 

where ©^(P^) = [P^ + m^j-i and Vpi^P"^) = [P^ + - impT{P'^)]-\ The 27r form factors 
of the physical vector mesons are 

■^loit-k{,P ) F'ujj-Kiri.P ) ~l~ ^{.P ) F^pimri^f ) 

Fp.AP') = -7=^^ {FpMP') - e(P')P.,..(P')) . (52) 





The form factors in the isospin basis have been described in Sec. 0. The dependence of the 
electromagnetic coupling strengths fp^ and /^^ on the mixing amplitude must be accounted 
for in a fit to data. (See, for example Ref. |^.) Here we limit our consideration to aspects 
of the 27r decay. Since the mixing amplitude e is part of a small second order effect in 
FpTrn{P^), it is the resonant u contribution, and in particular the physical coupling constant 
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g^^^^ = Fi^T^Tr{—m'^), which carries the leading information on the mixing amphtude. From 
Eq. (^2|) we may write 

by making use of the fact that u and uj have the same mass at the present level of approx- 
imation. Given that FpjT^T,{—m'^) is well determined by the known gpT^n, knowledge of gujjmr 
is necessary before e can be deduced from a value of g^j-Kn extracted from data analysis. 

For the calculated total mixing amplitude at the u mass shell, we find 11'"^ (—m^) = 
— (2.6 + 2.3i) X 10~^ GeV^. The real part is substantially below the experimental value 
of -(4.52±0.6) X 10-3 GeV^ from analysis of the e+e" — > TT+TT data. This latter value relies 
upon assumptions in the analysis that lead to the cancellation of the direct ujj — > tttt 



mechanism with Im 11'"^. A recent re-analysis has pointed out that the corrections to 
these assumptions lead to an increase of about a factor of three in the uncertainty of the 
extracted real mixing amplitude in the absence of theoretical limits on g^jj-w-K- Our result for 
Re 11''^ (—m^) is just within such an allowed range. The pion loop contribution to the p self- 
energy is found to be n''''(— m^) = —(0.31 + l.li) x 10"^ GeV^ giving a mixing parameter 



-m 



2 ' 



0.0303, arg(e) = 147.2° from Eq. (|T|). 



The influence of the direct coupling term g^jj-wK within the present model can be assessed 
as follows. With Eq. (|4l|) for e, Eq. (|53D becomes 

^ _Ren^ g^.^^WP - ilm IjP^gp,^^ 

where we write FpjT,T,{—m'^) ^ gpj-Kn since the shifted mass-shell point is a minor correction 
to this quantity. Although the amplitudes n^'^(P^) and Wp{P'^) here are to be evaluated at 
— m^, the arguments below are not altered significantly by changes in the mass-shell point 
compatible with the p — uj mass difference. The second term above is usually eliminated 



by making two assumptions Firstly, Rell^^ is taken to be zero. Apart from small 

adjustments due to the momentum dependence, this ignores the p — u mass difference given 
in Eq. (|^). Secondly, Im 11'"^ is assumed to be Im 11''^ scaled down by gui!TT-K I gpj-Kir- This 
would be exact in the present model if the pvrTr and utttt form factors were characterized by 
the same momentum dependent shape. As we have seen from Figs. ^ and |^, this is only true 
in a qualitative sense. As a result of the relevant loop integrals for the self-energies, there 
is an approximate scaling behavior evident in Figs. ^ and |^. Quantitatively, we find this 
second assumption to be in error by about 17% at the mass shell. To assess the net effect 
of these approximations, Eq. (|54D may be expressed as 



_ (Re n^-) C 



where the denominator is determined by the p mass and width via Eq. (|46|). The departure 
of the correction factor C from its conventional value of unity is a measure of the accuracy 
of the approximations mentioned above. Within the present model we find \C\ = 0.78 and 
arg(C) = 25°. A significant part of this correction is generated by (/o^jTrvrRe and thus is 
a measure of the competing influence of the usually neglected direct ujj tttt process. The 
trend is for a deduced value of Re II'"^ to be some 20% larger than what would otherwise 
be needed for the same gunrn- 
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VII. DISCUSSION 



In this work, we have found that a number of properties of the p — uj system can be 
described by a phenomenological treatment of dressed quark degrees of freedom as imple- 
mented by an effective field theory model applied up to one pion loop. A distinguishing 
feature of the present approach is that the finite size qq substructure of the mesons is in- 
cluded. This substructure for the pion is correlated with the momentum-dependent scalar 
self-energy dressing of the quark through chiral symmetry considerations. The resulting 
dynamical nature of the mass function in the pion propagator is included in the calculated 
pion loop contributions to both the mixed — u self-energy and the diagonal self-energy 
in the p channel. The latter produces a reasonable value for the p — uj mass difference if 
the influence of the qq substructure effects on the pion propagator are included. If we use 
the point particle limit of the pion propagator, this mass difference becomes negligible. The 
total p^ — UJ mixing here consists of a quark loop and a pion loop part and the latter is found 
to contribute 20% at the mass-shell. For spacelike momenta appropriate to mixed p — uj 
exchange for the NN force, the pion loop mixing contribution is not significant and the total 
mixing amplitude remains much too small to explain the empirical CSB NN force. 

The calculated values of QpTm and are used to fix two parameters that we have 
introduced to minimally add the vector mesons to a previously developed parameterization 
of pion physics that is the basis of the quark- level approach we employ. The isospin breaking 
quantities calculated here, i.e. the p^ — uj mixing amplitude and the coupling constant Quh-k-k) 
scale in an approximately linear way with the quark bare mass difference and we have 
presented our results for the typical value Sm = —4 MeV. 

We do not consider the contribution of p° — mixing to the meson masses because 
it is an effect of second order in the mixing amplitude. For a similar reason we do not 
consider the one pion loop correction to m^. In our pion loop calculations, we have used 
the empirical value of Qpt^t^ to set the strength of the p'jinr vertex even though the empirical 
coupling constant refers to the physical p° which is a linear combination of ^ind uji. For 
a calculation that is strictly of lowest order in the mixing amplitude, one should use the 
model results for the mixing to deduce a value for gpjnn and iterate the calculation until 
there is self-consistency. However consideration of Eq. (p2D shows that the correction is a 
small second order effect and we do not make it. 

In this work the coupling constant QujjTTn is found to be gujjmT = 0.105. We examine the 
correction that this direct process brings to the relation between the u; — >■ tttt decay and 
p^ — UJ mixing amplitude. The trend is for the necessary real mixing amplitude to be some 
20% larger than what would otherwise be the case. There is also a phase consequence and 
specific details can be found in Eq. (^). 

In a future work some of the approximations we make could be removed and some 
reduction could be sought in the degree of phenomenology employed. Although the strength 
of the single scalar function we use for the vector meson BS amplitude is determined by the 
mass-shell normalization internal to the model, the range of the gaussian form is described 
by a free parameter constrained by the calculated QpTr-K- We believe the result to be physically 
sensible because, in this same general framework, it leads to the prediction gf^^p = 0.50 
which is to be compared to g^^p = 0.504. Nevertheless, it is of interest to check these results 
with a vector meson amplitude that solves a BS equation and allows more than just the 
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canonical Dirac matrix covariant. Such numerical solutions are available [|26|. However it 
is difficult to accommodate the consequences of the complex values of quark loop momenta 
produced by mass-shell values of meson momenta in the Euclidean metric of this model. An 
ansatz for treatment of the ladder BS equation developed in recent works p7| , |37| goes a long 
way to overcoming this obstacle. Those works provide conveniently expressed BS amplitudes 
with parameters set by pion and kaon properties. One such approach allows a number of 
independent covariants selected by the dynamics and the predicted amplitudes for p and uj 
do conffim that the single canonical covariant used in the present work is strongly dominant. 
The other approach also provides a convenient form of the bare mass dependence of 
parameterized quark propagators that is somewhat more general than what we consider here 
and it would be of interest to explore the consequences. 
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APPENDIX: MESONS FROM QUARKS 

Hadronization of the quark action of Eq. (|1]) can be implemented through a functional 
change of field variables in the path integral for the generating functional. The meson sector 
follows directly from a standard Fierz reordering of the current-current term of the action 
in Eq. (|^) to obtain 

lf^[x)D{x-y)j;{y) = -\j\x,y)D{x~y)j\y,x), (Al) 

in terms of bilocal currents J^(x, y) = q{x)A^q{y). The A^ are direct products of Dirac spin, 
SU{2) flavor, and SU{3) color matrices, viz. 

A» = i ^OV. ^0V7.) S ^) » [lie, ^A") . (A2) 

We are here interested in the color singlet qq sector. There is however, a color octet qq sector 
in Eq. (|A^). An alternate color Fierz reordering may be used to eliminate the color octet qq 
sector in favor of color anti-triplet qq and qq pairings which find a natural role as constituents 
of baryons . However only the color singlet meson modes from Eq. (^^) are of concern 
here. The part of the integrand of the partition function Z = N J DqDq exp{—S[q, q]) that 
arises from the right side of Eq. (^1|) may be expressed as the integral 



exp 



^{J,DJ)^=N'J DB exp\^-^{B,D-'B)-iJ,B) 



(A3) 



where the bracket notation denotes the relevant space-time integration and summation 
over internal degrees of freedom. The net result is the reformulation Z =N" J DBDqDq 
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exp{—S[q,q,B]) where the new action has no four fermion term but instead contains 
(J, B)= J d'^x,y q{x)A^B^ {x,y)q{y), a hnear couphng of an appropriate set of bilocal auxil- 
iary fields B^ to quarks. To preserve the symmetries of the original action, it is necessary 
that the fields B^{x,y) have the same transformation properties as the currents q{y)A^q{x). 
These auxiliary fields are then bosons that take on the dynamics of qq correlations. For 
example, the u meson of the model will be dominated by the component associated with 
oc ^7^. With the action now quadratic in the quark fields, they may be integrated out by 
standard methods. The partition function then has the representation Z = J\f f DB e^"^''^' 
in terms of a bosonized action given by 

S[B] = -TrLn f ^ + m + A'B'] + f d'xd'y ^'(^' ^)^'(^' ^) . (A4) 
L J J 2D[x — y) 

An effective action for the propagating meson modes is defined by expansion in 
fluctuations about the saddle point of the action (which is equivalent to the classi- 
cal vacuum). From SS[Bo]/SBo = 0, these classical configurations given by Bq^t) = 
Z}(r)tr[A^(^ + m + A^i3^(r))~^] are related to nonlocal vacuum condensates and provide 
a quark self-energy S(r) = A^B^^r). This is equivalent to the ladder Dyson-Schwinger equa- 
tion in the bare vertex or rainbow approximation 

= J J2^^''^P - ^^y^^ .7-g + m + E(g) T^^- 

iFiom the general form =i ^[yl(p^) — 1] +B{p'^), the dressed quarks now have a 
dynamically-generated mass function M{p'^) = {B{p'^) + ■m)/A{p'^) and become the con- 
stituents of the meson modes. 

The meson modes are defined by expansion of the action about the saddle point config- 
uration of each auxiliary field. That is, one sets 

A^B%x, y) = S(x -y)+ i^i^f ■ vf (x, y) + ii^uj^{x, H , (A6) 

where 7f(x, y) + ■ ■ ■ are the new field variables, to produce the meson action from 

S[k,P,uj---] = S[B]-S[B^]. (A7) 

This will consist of terms at the quadratic and higher order level in the meson fields. Contact 



with effective local field variables may be made in the following way |Q . For each fiuctuation 
field 0, the free inverse propagator is given by = 5'^S/5(jP' at = 0. The free equation 
of motion V~^V = produces the ladder Bethe-Salpeter equation for the amplitude F. 
The set of eigenfunctions produced by V'^Tn = an^n may be used to expand the bilocal 
fluctuation fields thereby discretizing the continuous internal degree of freedom in terms of 
known amplitudes, leaving a local field as the remaining variable. Thus, for example, one 
arrives at the representation 7r(g; P) = J2n ^niQ] P)'^n{P) where q is the internal momentum 
conjugate to x — y, and P, conjugate to (x + y)/2, is the total momentum of the meson 
mode described by the local variable 7r„. With truncation to the lowest mass mode in each 
case, the result is Eq. (13). 
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FIG. 1. The quark-loop mechanism for puj mixing. 
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FIG. 2. The quark loop mechanism for pwir couphng. 
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FIG. 3. Dependence of the pyrvr form factor upon the p momentum with both pions on the 
mass-shell. The dot represents the p mass-shell. 
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FIG. 6. The real and imaginary parts of the pion loop contribution to the p self-energy. 
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FIG. 7. A graphical determination of the p — u) mass difference and the influence of the qq 
composite nature of the tt and p propagators as described in the text. 
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FIG. 8. The real and imaginary parts of the pion loop contribution to the — ui mixing 
amplitude. 



28 





FIG. 9. Various qq substructure effects on the pion loop contribution to p — oj mixing as 
described in the text. 
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FIG. 10. The quark loop (QL) and pion loop (PL) contributions to the — uj mixing amplitude. 
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